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Abstract 

Raman  scattering  studies  of  a  variety  of  MBE-grown  (GaAl)As  superlattices 
are  presented.  Folded  acoustic  phonons  appear  as  doublets  in  the  Raman 
spectra.  Their  frequencies  are  accurately  predicted  by  several  models, 
including  an  approximate  solution  of  an  elastic  continuum  model  through  a 
perturbation  approach.  Scattering  intensities  of  the  folded  acoustic  modes 
are  predicted  by  a  photoelastic  continuum  model.  Calculations  on  a  layered 
dielectric  continuum  provide  information  about  anisotropy  of  optical 
phonons.  Linear  chain  model  calculations  indicate  that  optical  phonons  in 
binary  superlattices  are  largely  confined  to  alternate  layers.  Peaks  in  the 
Raman  data  are  Identified  with  the  resulting  quantized  optic  modes.  It  is 
shown  that  Raman  scattering  has  the  potential  to  provide  structural 
information  similar  to  that  which  can  be  obtained  by  x-ray  diffraction. 
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I.  INTRODUCTION 


Semiconductor  superlattices  (SL's)  have  been  of  interest  for  nearly 
fifteen  years  now,  an  interest  stimulated  by  the  development  of  techniques, 
such  as  MBE  and  MOCVD,  capable  of  growing  crystals  of  alternating  layers  of 
two  semiconductors  with  nanometer  spatial  periods.  Potential  device  applica¬ 
tions  have  spurred  studies  of  the  electronic  properties  of  doped  and  undoped 
SL's  and  quantum  well  heterostructures.  Much  less  attention  has  been  paid  to 
the  vibrational  properties,  either  theoretical  or  experimental.  Calculations 
have  been  made  for  transverse  phonon-polariton  branches  in  a  GaAs-AlAs  SL* ) , 
and  discrete^)  and  continuum^)  models  have  been  used  to  describe  phonon 
folding  and  to  estimate  surface  mode  behavior.^ )  The  creation  of  a  gap  in  the 
acoustic  phonon  dispersion  by  the  SL  was  first  confirmed  by  acoustic  transmis¬ 
sion  experiments.5)  Raman  scattering,  a  natural  probe  for  SL  phonon  studies, 
has  been  the  most  frequently  employed  experimental  tool. 

Manuel  et  al.'’ '  studied  the  resonant  Raman  behavior  of  LO-phonons  near  the 

energy  gap  of  four  GaAs-A^Gaj^As  samples  with  63A  <  d  <  20QA  and  x  ~  0.25. 

They  also  derived  a  two-dimensional  result  for  the  resonance  behavior.  An 

extensive  study  was  made  by  Barker  et  al.^)  of  Raman  scattering  and  infrared 

absorption  in  a  series  of  very  thin  layer  GaAs-AlAs  samples.  They  calculated 

the  effects  of  phonon  folding,  and  made  comparisons  between  alloy  samples  and 

7  8) 

a  number  of  thin  layer  SL's#  Sai-Halasz  et  al»  ’  '  noticed  new  peaks  appearing 


i 

Is 


between  the  LO®8^8  and  LoAlGaAs  modes  in  resonant  Raman  spectra  of  GaAs- 
Al'25Ga.75A8  samples  which  they  attributed  to  umklapp  processes.  These  new 
peaks  were  shown  to  be  similar  to  E(L0)  phonons  by  Merlin  et  al.®)  The  first 
clear  phonon-folding  in  light  scattering  experiments  was  seen  by  Colvard 
et  al.5)  Later  Raman  studies  report  scattering  from  both  acoustic*®’**)  and 
optic*®)  phonons,  including  resonant  enhancements  at  quantized  exclton  levels.*^) 
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Brillouin  scattering  experiments  have  been  performed  by  Sapriel  et  al.*^ 
Indications  of  optical  phonon  folding  have  been  reported  by  Jusserand  et  al.*^ 

We  present  here  a  comparison  of  experimental  Raman  scattering  results 
from  a  variety  of  SL's  with  several  theoretical  models,  including  a  pertur¬ 
bation  calculation  which  reproduces  the  essentials  of  the  acoustic  phonon 
spectrum,  such  as  sound  velocity  and  gap  widths.  We  show  that  Raman  data  can 
provide  information  on  the  layering  parameters  heretofore  gathered  primarily 
from  x-ray  diffraction.  A  model  based  on  bulk  photoelastic  behavior  is  shown 
to  predict  the  Intensity  of  the  observed  folded  acoustic  phonon  peaks.  The 
optical  phonons  are  seen  to  be  best  described  by  quantization,  as  in  the 
electron  case,  rather  than  by  zone  folding. 

In  the  next  section,  several  models  of  SL  phonon  behavior  are  outlined. 
Section  III  presents  a  calculation  of  light  scattering  intensities,  relating 
the  folded  phonon  and  Brillouin  scattering  intensity  to  the  bulk  photoelastic 
constants.  Sections  IV  and  V  present  the  experimental  results  and  discussion. 


II.  PHONON  MODELS 


I 

Several  different  approaches  to  modelling  superlattice  phonons  have  been 
found  to  be  useful,  each  having  certain  advantages.  A  linear  chain  model 
embodies  the  discrete  nature  of  the  crystal  structure  and  provides  information 
about  evanescent  behavior  of  optical  phonons.  It  is  most  useful  in  frequency 
regions  in  which  the  bulk  phonon  dispersion  is  non-linear.  At  low  frequencies 
and  long  wavelengths  the  lattice  can  be  treated  as  an  elastic  continuum,  in 
which  case  a  Kronig-Penney  type  of  approach  yields  a  convenient  analytic 
expression  for  phonon  folding.  The  layer  thicknesses  in  this  case  are  not 
restricted  to  integral  multiples  of  the  monolayer  spacing.  Alternatively,  at 
these  same  frequencies,  the  superlattice  can  be  considered  as  a  perturbation 
on  an  average  elastic  background.  This  approximate  approach  yields  explicitly 
the  gap  widths  and  sound  velocity  for  an  arbitrary  composition  profile. 

These  models  will  be  outlined  below,  in  addition  to  a  dielectric  continu¬ 
um  model  which  allows  an  estimate  of  the  effects  of  long-wavelength  electric 
fields.  Of  course,  one  wants  a  comprehensive  theory  of  superlattice  effects 
which  correctly  deals  with  all  microscopic  short  and  long  range  forces.  Such 
a  theory  is  especially  needed  in  dealing  with  the  spectrum  of  optical 
phonons.  Calculations  using  such  a  theory  have  recently  been  made  by  Yip  and 
Chang. 

A.  Linear  Chain  Model 

This  model  has  been  applied  to  binary-binary  superlattices  by  Barker  et 
al.^)  it  is  valid  in  bulk  crystals  for  phonons  propagating  along  [001],  where 
planes  of  atoms  move  as  a  whole,  and  the  longitudinal  and  transverse  vibm- 
tions  are  decoupled.  For  longitudinal  modes  a  single  nearest-neighbor  spring 
constant  provides  the  only  free  parameter,  whereas  two  different  nearest- 
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neighbor  spring  constants  provide  the  anisotropy  necessary  to  describe  trans¬ 
verse  vibrations. ^  These  force  constants  are  allowed  to  be  different  in  GaAs 
and  AlAs  layers,  and  are  fit  to  bulk  values  of  either  the  L0(T)  phonon  frequency 
(or  T0(r)  and  TA(X)  for  transverse  modes)  or  the  sound  velocity  viong  (or 
vtrans  an<*  ^®(r)  frequency).  To  fit  both  optic  and  acoustic  regions,  more 
distant  neighbor  interactions  must  be  included. 

By  utilizing  bulk  phonon  solutions  (propagating  or  evanescent)  within  the 
layers  and  matching  these  at  the  interfaces,  arbitrary  layer  thicknesses  can 
be  considered  by  solving  only  a  4x4  determinant.  The  longitudinal  case  will 
be  described  here;  the  extension  to  anisotropic  force  constants  for  the 
transverse  case  is  straightforward.  Figure  1  shows  the  model  used.  Within 
each  layer 


inline]  -  -  K  (21^  [me]  -  U2[me]  -  U2[(m-l)eJ)  (1) 

M2U2 [me ]  -  -  K  (2U2[meJ  -  UxKm  +  l)ej  -  UjJmeJ)  (2) 

* 

where  U(me)  is  the  displacement  along  z  of  the  atom  at  z  ■  me,  and  e  is  the 
monolayer  spacing  e  -  a/2.  Taking  ±a  as  the  phonon  wavevector  qz  in  GaAs, 
this  yields 

(M  w2-  2K)  (M?u>2-  2K)  -  2K2 

cos  ae  - - - -  .  (3) 

2K 


k'/A 


PW 
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V 


fcV 


Defining  Y±  =  l^/Uj  for  ±a,  we  find 


-k  a  +  .llae) 

t  2 

M2«  -  2K 


(4) 


The  above  equations  are  for  GaAs;  similar  expressions  describe  AlAs  with 

wavevector  qz  *  ±8  and  6^5  U^/Uj.  ^q.  O)  *s  Plotte<*  *n  Pig*  2,  which  shows 

both  the  real  and  the  imaginary  part  of  qz*  The  imaginary  part  indicates  the 

rate  of  attenuation  of  a  vibration  at  a  particular  frequency.  For  y  -  2.83A,  a 

penetration  depth  of  one  monolayer  for  the  amplitude  U(z)  corresponds  to 

Im(qz)  -  0.32  */(y).  The  plot  indicates  that  at  GaAs  optical  mode  frequencies 

the  vibrations  extend  less  than  one  monolayer  into  the  AlAs,  where  they  are 

damped  optical  modes  with  Re(q)  -  2* /a.  At  AlAs  optical  frequencies  they  are 

even  more  localized  in  the  AlAs  layer,  decaying  into  the  GaAs  as  damped 

acoustical  modes  with  Re(q)  -  0.  The  result  of  this  confinement  of  the 

optical  vibrations  is  apparent  in  Fig.  3,  where  the  folded  superlattice 

optical  modes  are  seen  to  be  quite  flat. 

In  a  superlattice  the  displacements  become  Uj(z)  «  Aeiaz  +  Be-iaz  and 
iBz  -iSz 

Uj(z)  ■  Ce  +  De  .  At  an  Interface  the  forces  must  be  matched,  so  that 

each  wave  behaves  as  if  it  were  in  an  infinite  medium.  Thus  at  z  =  0: 

K[U2(0)  -  U^O)]  -  K*[U£(0)  -  UJ  (0)] 

and 

KfUjCe)  -  U2(0)J  -  K* [ UJ  (e)  -  U£(0)J. 

With  the  help  of  Eq.  (1)  this  gives 

U2(0)  -  4(0)  (7) 

and 

KlU^e)  -  4(0)]  -  K'{U[(e)  -  UJ(0)J  (8) 


(5) 

(6) 


which  are  the  discrete  analogs  of  the  requirements  that  displacement  and 


stress  be  continuous  at  an  acoustical  boundary.  Similarly,  at  z  -  ne  we  have 


U2(ne)  -  U£  (ne) 


KlU^ne  +  e)  -  U^neJJ  -  K*[U|(ne  +  e)  -  U|(ne)J. 


Imposing  the  periodic  requirement  that  UgL(z+d)  -  UgL(z)  eiqd,  where  q  is 
the  super lattice  wavector,  and  setting  layer  thicknesses  d^  +  d2  “  d,  gives 
the  set  of  equations : 


Y+A  +  Y_B  -  «+C  +  <S_D 


G+A  +  G_B  -  H+C  +  H_D 


-iad  iqd  iad  iqd  iSd  -iPd 

Y^Ae  e  +  Y_Be  e  =*  6+Ge  +  6_De 


-iad.  iqd  iad  iqd  iBd  -iSd 

G+Ae  e  +  G_Be  e  -  H+Ce  L  +  H_D  e  1 


(11a) 


(lib) 


(11c) 


(lid) 


where 


G±  i  K(e±icte-1) 


H+  s  K'(e±iBe-l) 


(12a) 


(12b) 


The  solution  of  these  equations  is  given  by 


7 


cos  qd  -  cos  ad^  cos  0d2 


2(Y+Y_H+H_+  «+«_G+G_) 

(Y_6_g+H++  Y+6+G_H_-  Y_6+G+H_-  Y+«_G_H+)  ®in  °dl  8ln  Pd2 


(13) 


In  the  transverse  case,  an  additional  term  appears  in  the  sine  multiplicand. 
This  solution  is  shown  in  Fig.  3  for  a  (5,4)  superlattice.  The  phonon 
branches  above  the  first  but  below  about  200  cm-*,  although  strictly  optical 
branches  of  the  SL,  are  commonly  referred  to  as  folded  acoustic  branches  due 
to  their  origin  in  the  bulk  acoustic  phonons. 


B.  Elastic  Continuum  Models 

The  phonon  dispersion  is  approximately  linear  below  about  100cm-*  for  LA 
phonons  and  below  50  cm-*  for  TA  phonons,  and  in  this  region  the  crystal  can 
be  treated  as  an  elastic  continuum.  Vibrations  in  layered  elastic  media  have 
been  considered  by  Rytov*^,  and  his  results  have  been  shown  to  apply  to 
acoustic  phonons  in  superlattices.^  The  dispersion  is  given  by 

“di  “d?  1  +  x2  wdt  “d? 

cos  qd  -  cos(- — )  cos(— -~-) - j- —  sin  (— — -)  sin  (— — )  (14) 

V1  V2  Z  V1  2 
which  is  a  solution  of  the  wave  equation 

(15) 

Here  the  sound  velocity  is  v  *  (K/p)*^  in  each  layer,  and  <  =  (>2V2^  lvl* 

In  the  limit  of  large  wavelengths,  the  superlattice  sound  velocity  v<^  -  <*»/q 


becomes 


If  Che  elastic  constants  are  equal  in  the  two  layers,  this  reduces  to 


VSL  “  ^<v_1>rms^-1 *  U7> 

The  above  Rytov  solution  is  valid  for  a  square-wave-like  compositional 
modulation  with  sharp  interfaces.  Its  results  can  be  recovered,  and  extended 
to  arbitrary  composition  profiles,  by  a  perturbation  approach.  We  replace 
Eq.  (15)  by 


[(K  +  AK)  j^]  +  (Pq+  Ap)  w2U  =  0 


(18) 


where  KQ  and  pQ  are  the  zeroth  Fourier  components  of  the  elastic  constant  and 
density,  and  AK(z)  and  Ap(z)  describe  the  compositional  variation.  They 
depend  upon  the  local  A1  concentration  x.  The  lowest  order  equation  is  then 


K 

P 


o  afu 
o  3z2 


+ 


ui2U 


0 


(19) 


o 

which  gives  a  virtual  crystal  sound  velocity  v*  =  ^o7po*  ’^1^8  agrees  with 
Eq.  (17)  for  equal  elastic  constants. 

We  can  make  the  expansions 


U(z)  -  I  U  eiqz 


(20a) 
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Ap  =  l  Pce 

S  * 


(20c) 


where  S  is  a  reciprocal  lattice  vector,  S  -  2xm/d,  and  q  ranges  over  the  bulk 
Brillouin  zone.  Eq.  (18)  then  becomes 


(voV-  „2)  Uq+  l  A-  l,(,-S)  Ks-  .2pg]  Uq_g  -  0. 

S  o 


The  ur  in  the  perturbation  term  in  Eq.  (21)  can  be  approximated  by 


u2-  V2  | q(q— S) | 


This  allows  the  perturbation  to  be  written  as  Vg  q(q-S),  where  we  define  Vg  £ 

_  t 

(Kg  +  V  p  )/p  .  The  -(+)  sign  applies  if  q(q-S)  is  positive  (negative). 

^  O  O  o 

We  keep  only  the  terms  in  the  sum  which  correspond  to  states  that  are  strongly 
mixed  by  the  perturbation,  i.e.,  for  which 


lvoq2“  vo(q  “  S)2I  *  |q(q-s)vgl 


This  is  true  for  the  states  and  Uq_g  where  S  *  2q.  Considering  only  these 
two  terms,  we  can  simplify  Eq.  (21)  and  obtain 


»2-~  [q2+  (q-S)2]  ±  [-1  (q2-  (q  -  s)2)  +  |V.|2q2(q-S)2]1/2, 


which  for  exact  degeneracy,  S  ■  2q,  becomes 


«2-  vo2q2±  |VS|  q  (q-S). 


'.V.V.Y -.-.V. 


If  we  set  q  *  -j—  then 


“2  ■  (vo2±  lVsl>  ^ 


gives  the  frequencies  on  either  side  of  a  gap,  at  the  SL  zone  edge  for  m  odd 
and  zone  center  for  m  even.  Taking  pg  for  a  square  wave,  and  assuming  Kg  « 
(i.e.  K  equal  in  the  two  layers),  we  find 


v  (p.  -  P  )  nmd, 

| v_ |  -  |-2— 2 - —  sin  (__JL)i 

'S'  '  p  it  m  d  1 

o 


where  pa  and  p^  are  the  densities  in  the  two  layers,  and  where 


P *  d  ^dip„+  d„P. ). 

o  1  a  L  b 


The  expected  gaps  are  found  by  noting  that 


lVs'  ilmj 

*  d 
o 


which  gives,  again  for  constant  K  and  sharp  boundaries, 


Dielectric  Continuum  Model 


For  layers  thick  enough  that  a  local  dielectric  constant  is  a  well 
defined  quantity,  the  optical  phonon  behavior  can  be  surmised  by  treating  the 
SL  as  a  layered  dielectric  continuum.  The  dielectric  constant  can  be 


represented  in  the  binary  layers  by 


e(u>) 


e.(l  + 


2  2 

us^r  ui  -  io>r 


(31) 


which  for  zero  damping  becomes 


2  2 
to  —  at 

e<“)  "  e-  (“J - * 

0)1-  0) 

T 


In  the  alloy  we  can  use  the  factorized  form^^ 


e(u>) 


(o)^-  o >2)  (ui^2—  “2> 

<»„-  <*>2)  “2) 


(32) 


(33) 


where  u>L  and  u>T  are  the  long  wavelength  longitudinal  and  transverse  phonon 
frequencies.  The  subscripts  1  and  2  refer  to  the  composition  dependent  two¬ 
mode  behavior  of  the  alloy.  Maxwell's  equations  can  then  be  solved  for  polari 
zation  waves  propagating  parallel  or  perpendicular  to  the  layers  in  the  q  ♦  0 
limit.  This  has  been  done  in  the  absence  of  phonon  dispersion  by  Rytov.*®^ 

The  results  can  most  conveniently  be  summarized  by  noting  that,  for  long 
waves,  the  superlattice  behaves  like  a  homogeneous,  anisotropic  crystal.  The 
average  dielectric  constants  are  found  to  be 


£*,y-  J  •  d'1(dlV  d2V  (34> 


v  t*'1  r1-  dE1E2(diV  d2eir‘  (35) 


with  and  62  given  by  Eqs.  (32)  or  (33) 


The  zeroes  and  poles  of  these  dielectric  constants  occur  at  the  longi¬ 
tudinal  and  transverse  frequencies  of  the  composite  layered  structure, 

corresponding  respectively  to  the  vanishing  of  the  average  displacement  D  or 

-  ♦ 
electric  field  E.  For  polarization  in  the  plane  of  the  layers  and  q  along 

A 

the  z  direction,  the  poles  of  Eq.  (34)  give  TO  modes  of  E-symmetry  at  the  bulk 
TO  frequencies.  For  polarization  parallel  to  the  superlattice  axis,  the  zeroes 
of  Eq.  (35)  give  LO  modes  of  or  A^-symmetry  at  the  bulk  LO  frequencies. 

In  both  cases  the  polarization,  and  electric  field  where  non-zero,  exists  only 
in  the  layer  at  whose  normal  frequency  the  vibration  occurs.  If  phonon 
dispersion  is  considered,  these  modes  are  the  quantized  optic  phonons  due  to 
SL  zone  folding. 

For  3  in  the  layer  plane  two  new  modes  appear,  a  ^(TtO-like  mode  at  the 
poles  of  Eq.  (35)  and  an  E(LO)-like  mode  at  the  zeroes  of  Eq.  (34).  The 
solutions  are  given  by 

(36a) 

(36b) 

with  Eqs.  (32)  or  (33)  for  and  C2>  These  frequencies  occur  between  the 
bulk  and  u>T.  For  qd  >  0  the  polarization  decays  exponentially  away  from 
interfaces,  but  in  the  thin  layer  limit  it  is  nearly  constant  throughout  the 
superlattice.  Figure  4  plots  the  frequencies  of  these  modes  versus  layer 
thickness  for  a  GaAs-AlAs  superlattice.  Note  that  if  d^  »  d£,  the  LO-mode 
can  occur  at  the  TO-frequency  and  vice-versa. 

Solutions  will  also  exist  at  the  bulk  LO  and  TO  frequencies  for  modes 

A 

propagating  along  the  layers  but  spatially  quantized  along  z.  Ignoring 
polariton  effects,  these  modes  will  be  similar  to  the  quantized  bulk  phonons 
mentioned  above  at  the  wavevectors  accessible  to  Raman  scattering. 


III.  Photoelastic  Mechanism  for  Light  Scattering 

At  laser  energies  far  from  resonance  with  superlattice  electronic  transi¬ 
tions,  light  scattering  from  the  folded  acoustic  phonons  can  be  treated  as  a 
coherent  sum  of  scattering  within  bulk-like  layers  due  to  the  photoelastic 
effect.  The  photoelastic  coefficient  is  taken  to  be  a  function  of  position  z 
its  value  depends  on  the  local  A1  concentration  x: 

P(z)~Sx(z)*  (37) 

If  the  interfaces  are  abrupt,  this  will  be  a  square  wave: 

P  0  <  z  <  d. 

a  —  l 

P(z)  -  (38) 

Pb  dj<  z  <  d  -  d1+  d2 

The  strain  due  to  an  LA  phonon  will  cause  a  proportionate  change  in  the 
susceptibility : 


«X(Z)  -  P(z)  .  (39) 

f*  h 

Light  will  be  scattered  by  the  component  of  this  susceptibility 
fluctuation  that  conserves  momentum,  q  -  £  -  k  ,  as  in  ordinary  Brillouin 
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scattering 


I(o))  «  <  |6xq|2  V 

This  component  is  given  by 


(40) 


V  C  e-'^PU)  ^  d. 
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In  the  Rytov  model,  away  from  gaps  in  the  phonon  dispersion  curves,  this 
integrand  is 


<*X(z) 
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0  <  z  <  d. 
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where  a  =  — ,  0  »  —  and  v,  2  are  bulk  sound  velocities.  For  simplicity  the 
1  v2 

difference  between  a  and  0  can  be  ignored,  a  “  0  *  q  and  U._  *  U  .  .  P(z)  can 
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be  separately  expanded 


P(z)  -l  P  e  "  ,  S  =  ~r- 

L  m  ’  m  d 


allowing  the  Raman  tensor  to  be  written 


fiX  *  ^  Pm  “  SiJ  Un  <5  * 
qm  m  q— S 


The  ±  m  terms  in  this  sum  describe  Raman  scattering  from  folded  LA 
phonons  at  «,  ■  |q  +  Sm|vSL  with  intensity 


I  «  ut  (n  +  1)  |P  I  . 
m  m  m  '  nr 


(Here  is  the  Bose  factor).  For  m  *  0  this  gives  the  Brillouin  scattering 
cross  section  in  a  superlattice  at  u  ■  q  vci  •  For  the  square  wave  case  of 


Eq.  (38) 


Pn-  4  (d,  P  +  d„P,  ), 


WWv! 


r 


i 


* 


£ 


i(Ph"  Pa) 

V  - 2xm -  [1-  expC-iZumd^d)], 


m  *  0, 


giving 


.  (P.-  Pa)  nurd, 

|P  |  ■  -  sin  (— — )  m  *  0. 

1  m1  mu  '  d  ' 


For  propagation  along  z  and  backscattering  from  a  (001)  face,  only  LA  phonons 
can  be  seen,  these  through  the  component  pXX2Z  =*  p^-2  in  cubic  crystals.  Thus 
the  scattered  folded  phonon  intensity  can  be  compared  with  that  for  the 
Brillouin  scattering: 

(Pb2-  P«12)2  si*i2(mwd.  /d) 

V^Brill*  ~2  T1  n  (49) 

P0 


The  factor  n  accounts  for  the  frequency  difference  between  u>m  and 

n  -  tU)m(nm+  1)]/[(1,o^no+  For  wm«  kT,  n  •  l. 

Ren  and  Harrison  have  calculated  photoelastic  constants  for  GaAs  and 
A1As.2(^  As  defined  here,  their  values  give  P12  =  0.48  in  GaAs  and  0.05  in 
AlAs  for  the  change  in  x  with  strain.  For  dj  =  d2  =  d/2,  this  gives 


Il/IBrill  "  °*3  * 


IV.  Experimental 

All  samples  used  in  this  study  were  grown  by  molecular  beam  epitaxy  on 
[001]  oriented  GaAs  substrates.  Periods  range  from  20A  to  20QA,  and  SL 
thicknesses  vary  between  0.5y  and  4y.  Specific  SL's  will  be  referred  to  by 
their  composition.  The  notation  used  for  an  AlyGaj_yAs-AlxGaj_xAs  SL  will  be 
(dj »y>d2 :x) »  wtiere  dj  and  d2  are  layer  thicknesses.  If  no  value  of  y  is 


specified,  as  in  most  samples,  the  first  number  refers  to  pure  GaAs  layers 
(y  -  o). 

Light  from  a  cw  Ar+,  Kr+,  or  dye  laser  was  used  to  illuminate  the  samples 
which  were  held  in  a  He  cryostat  at  temperatures  between  2-300K.  Scattered 
light  was  analyzed  by  a  double  pass  3/4  meter  spectrometer  equipped  with  1800 
g/mm  holographic  gratings  and  was  detected  by  an  RCA  C31034A  photomultiplier 
tube  operated  in  phonon-counting  mode.  An  Apple  II-plus  computer  was  used  to 
collect  data  and  control  the  experiment.  Most  measurements  were  made  with 
light  incident  at  an  angle  close  to  Brewster's  angle,  polarized  in  the  plane 
of  incidence.  Scattered  light  was  collected  with  an  f/1.4,  85  mm  focal  length 
lens  whose  axis  was  normal  to  the  surface  of  the  sample.  Successive  rotations 
of  a  polarizer  in  the  scattered  beam  allowed  polarized  and  depolarized  spectra 
to  be  Interspersed  in  the  same  spectral  scan,  so  that  small  frequency  shifts 
between  the  two  spectra  could  be  detected. 

Pieces  of  one  (4lA,  8A:1)  sample  were  annealed  at  850°C  for  times  from 
0.25  to  16  hours.  During  the  anneals,  the  samples  were  placed  between  two 
pieces  of  GaAs  in  an  atmosphere  of  flowing  forming  gas  inside  a  quartz  furnace 
No  visible  surface  degradation  or  increased  surface  scattering  was  seen  after 
annealing  for  up  to  16  hours. 


Results  and  Discussion 


Acoustic  Region  (5-220cm  ) 


In  this  region  are  seen  the  most  striking  features  of  a  superlattice 
spectrum:  the  folded  acoustic  phonons.  They  appear  as  doublets  in  the 
spectrum  at  frequencies  from  5  to  100cm-1.  An  example  is  seen  in  Fig.  5, 
which  exhibits  three  folded  orders,  as  shown  in  the  inset.  The  data  are  taken 
from  a  (42A,  8A:.3)  sample.  Observation  of  small  peaks  close  to  the  laser 
line  is  greatly  aided  by  the  quality  of  the  MBE-grown  sample  surface.  Peaks 
as  low  as  7  cm-1  have  been  clearly  seen.  At  these  small  shifts,  the  intensity 
of  the  light  scattering  is  enhanced  by  the  thermal  factor  n  +  1,  which  at  room 
temperature  is  21.5  at  10cm-1,  compared  to  1.3  at  the  292cm-1  GaAs  LO  phonon 
frequency.  The  data  presented  here  include  this  factor,  which  must  be  kept  in 
mind  when  interpreting  relative  intensities. 

Figure  6  shows  some  further  examples  of  the  folded  doublets,  all  spectra 
taken  out  of  resonance  at  5145A.  The  widths  of  the  peaks  are  quite  small, 
usually  2  to  3  cm-1,  including  an  instrumental  linewidth  of  about  2cm-1.  The 
expected  contribution  to  the  width  due  to  absorption  broadening  of  q  is  less 
than  0.5cm-1.  The  acoustic  phonons  apparently  experience  a  well  defined 
average  period,  and  show  little  inhomogeneous  broadening.  Also  seen  in  the 
spectra  is  the  well  known  2TA  structure,  which  peaks  at  160cm-1  in  GaAs  and  at 
about  200cm-1  in  AlAs.  After  correcting  for  the  n  +  1  Bose  factor,  we  find 
that  the  intensity  of  the  first  folded  peak  varies  from  0.01  to  0.13  times  the 
allowed  GaAs  LO  intensity  in  the  samples  studied,  averaging  about  0.05.  Under 
resonant  enhancement  these  two  Intensities  can  become  comparable. 

As  pointed  out  in  the  previous  section,  the  folded  doublets  correspond  to 
scattering  from  acoustic  phonons  with  wavevectors  2trm/d  ±  q.  For  frequencies 
away  from  the  gaps,  i.e.,  where  the  acoustic  dispersion  is  linear,  these 
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longitudinal  modes  are  running  waves  with  equal  components  of  A^  and  B2 
character.  Light  scattering  from  the  Raman-active  Aj  component  gives  equal 
intensities  for  the  members  of  each  doublet. This  can  be  seen  in  the  lower 
curves  of  Fig.  6.  The  upper  curves  show  the  relative  intensities  of  the 
doublets  becoming  progressively  more  unequal  in  samples  with  smaller  periods. 
This  indicates  some  admixture  of  a  standing  wave  due  to  Bragg  reflection  and 
the  opening  of  a  gap  in  the  dispersion  curve.  Figure  7  illustrates  the  first 
zone-center  gap  in  a  (5,4)  superlattice.  The  solid  curve  is  from  a  calcula¬ 
tion  using  the  linear  chain  model,  which  fits  the  bulk  sound  velocity  of  GaAs 
with  the  spring  constant  being  the  only  adjustable  parameter.  The  data  points 
were  taken  from  a  sample  with  25 .5A  period  by  varying  the  laser  wavelength 
from  4579A  to  6764A.  In  such  thin  period  samples,  the  wavevector  transfer  q 
is  just  barely  in  the  region  where  the  dispersion  begins  to  flatten  out  into 
the  gap.  In  the  case  of  thicker  GaAs  layers,  i.e.  dj  >  d2,  the  upper  curve 
takes  on  B2  symmetry,  the  lower  one  A^  symmetry.  Because  in  acoustic  modes 
the  atoms  in  a  unit  cell  move  in  phase,  the  B2  mode  is  effectively  of  odd 
symmetry  under  inversion.  This  gives  a  vanishing  cross-section,  since  purely 
odd  modes  cannot  ordinarily  participate  in  Raman  scattering.  One  can  predict 
that  small  period  samples  having  d^  <  d2  will  display  the  opposite  asymmetry, 
the  lower  frequency  peak  having  mostly  B2  symmetry  and  thus  being  weaker. 

The  wavevector  transfer  in  backscattering  is  too  large  to  truly  probe  the 
folded  gap.  In  forward  scattering,  although  qz  can  be  made  small,  spurious 
laser  light  can  overwhelm  data  taken  at  small  frequency  shifts.  Another 
possible  way  to  study  the  gaps  is  to  make  samples  with  laige  enough  period  so 
that  qz  =■  4irn/XL  *  j  ,  allowing  scattering  from  the  zone-edge  gap.  This  would 
require  d  3  250A  at  XL  *  4579A  in  a  sample  with  n  3  4.5,  well  within  the  realm 
of  current  technology. 


The  data  presented  thus  far  are  all  on  longitudinal  phonons.  Some  folded 


transverse  phonons  have  been  seen  in  a  Brewster  angle  scattering  geometry. 

Figure  8  shows  an  example.  The  transverse  phonons  labeled  by  T  appear  in  both 
(x,x)  and  (x,y)  spectra.  They  have  E  symmetry  in  the  superlattice,  and  are 
seen  due  to  the  deviation  from  true  backscattering. 

A  comprehensive  confrontation  of  theory  vs.  experiment  for  the  folded 
phonon  frequencies  is  difficult,  considering  the  variety  of  periods,  layer 
thickness  ratios,  and  alloy  compositions  studied.  In  addition,  the  bulk  sound 

velocity  in  AlAs  is  not  well  known.  It  is  often  estimated  (as  here)  by 

2 

v  =  C/p,  using  the  same  elastic  constant  C  as  in  GaAs.  (The  values  of 
s  i 

elastic  constants  Cj  ^  =  12.5,  =  5.34,  and  -  5.42*  lO^dyn/cm^  are 

given  in  reference  20  for  AlAs,*  estimated  from  empirical  relations  and  typical 
III-V  values.  This  Cj  j  is  6%  larger  than  in  GaAs).  Figure  9  gives  a  plot  of 
the  folded  phonon  frequencies  observed  vs.  inverse  period.  The  shaded  bands 
indicate  the  frequencies  expected  for  the  first  three  foldings  of  an  average 
linear  acoustic  branch.  They  are  bounded  by  u>±-  v±(2xm/d),  where  v+  is  the 
sound  velocity  of  AlAs,  v_  that  of  GaAs.  The  midpoint  of  each  doublet  is 
expected  to  lie  within  such  a  band,  its  exact  position  depending  on  the  A1  to 
Ga  ratio  in  the  sample.  The  doublets  off  to  the  lower  right  represent 
transverse  phonons. 

The  doublet  splitting  is  seen  to  be  -  5cm- ^  in  all  observed  cases.  This 
is  consistent  with  the  wavevector  transfer  q  **  4nn/A^  for  backscattering. 

The  phonon  doublets  are  at  u>  =  Vg^  (Sm  ±  q)  giving 

Am  -  2vSLq  -  8irnvSL/AL 
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1C  is  felt  that  the  scatter  in  the  data  is  largely  due  to  uncertainty  in 
the  value  of  the  sample  period.  When  this  was  not  determined  by  x-ray 
diffraction,  it  was  estimated  from  growth  conditions*  Furthermore,  a  sample's 
period  can  vary  across  its  surface.  Figure  9  demonstrates  the  ability  of 
Raman  scattering  in  a  superlattice  to  provide  structural  information,  similar 
to  that  from  x-rays,  by  measuring  the  sample  period.  It  has  the  advantage  of 
providing  convenient  spatial  resolution  as  well. 

Another  feature  of  interest  in  Fig.  9  is  the  appearance  of  two  samples 
showing  m  -  1  and  m  *  3  phonon  doublets,  but  none  with  m  -  2.  Figure  10 
presents  data  from  sample  (4lA ,4lA  :0.3) .  The  inset  is  a  dispersion  calcu¬ 
lation,  using  the  Kytov  model,  with  the  observed  peaks  indicated.  Missing 
peaks  are  also  seen  in  the  data  of  Figure  11  on. a  (85A,88A  :0.3)  sample.  The 
inset  corresponds  to  the  lower  curve,  taken  out  of  resonance  at  XL  -  5145A. 
(Note,  incidentally,  that  the  data  are  taken  with  q  over  halfway  into  the 
superlattice  Brillouin  zone.)  These  data  lend  support  to  the  model  of  photo- 
elastic  scattering  presented  in  section  2.  Both  samples  have  nearly  equal 
layer  thicknesses;  Eq.  (48),  along  with  Eq.  (45),  shows  that  in  this  case  the 
scattering  vanishes  for  even  values  of  m.  The  x-ray  diffraction  data  in  Fig. 

12  were  taken  from  this  sample.  Here  also  the  second  satellite  is  missing, 
corroborating  these  results.  In  the  samples  studied,  the  ratio  of  the  inten¬ 
sity  of  the  m  ■  1  peak  to  that  of  the  m  =  2  peak  is  generally  consistent  with 
Eq.  45,  after  corrections  are  made  for  the  thermal  factor  in  Eq.  (45).  The 
intensities  of  the  second  and  third  peaks  in  Fig  5  fall  off  slightly  faster 
than  predicted,  presumably  due  to  the  sensitivity  of  sin^(wmdj/d)  to  the  exact 
layer  thicknesses,  especially  for  larger  m. 

To  see  the  effect  of  a  non-sharp  interface,  several  pieces  of  sample 
(4lA,8A:l),  which  contains  pure  AlAs  barrier  layers,  were  annealed  at  850°C  for 


several  intervals,  up  to  16  Hours.  This  causes  the  interfaces  to  diffuse.  1 
The  effect  on  the  folded  phonons  is  seen  in  Fig.  13.  The  intensity  decreases, 
while  the  frequency  and  width  (within  the  experimental  limits)  do  not  change. 
This  agrees  with  predictions  of  the  photoelastic  scattering  model.  The 
background  at  frequencies  below  that  of  the  phonons  also  decreases.  The  data 
were  taken  at  room  temperature  and  are  normalized  to  50  counts/s  at  265cm-* 
(the  TO(GaAs)  peak).  No  significant  increase  was  seen  in  disorder  activated 
scattering  (DALA  and  DATA). 

Under  resonant  conditions  the  intensity  of  the  folded  phonon  scattering 
can  be  significantly  altered.  An  example  is  the  upper  curve  of  Fig.  11,  taken 
near  the  Eq  +  Aq  gap  of  the  superlattice.  The  even-m  phonons  are  now  seen, 
but  the  peaks  are  asymmetric,  and  they  are  shifted  in  a  way  that  cannot  be 
explained  entirely  by  a  shift  in  q.  Resonance  of  the  laser  photon  with 
excitation  to  a  particular  SL  state,  which  is  an  appropriate  combination  of 
bulk  propagating  and  evanescent  states  in  the  two  layers,  invalidates  the 
assumptions  of  Section  III.  The  asymmetry  suggests  a  coupling  between  the 
discrete  phonons  and  a  background  continuum.  Although  not  pursued  here,  other 
examples  of  such  a  coupling  have  been  seen,  as  in  Fig.  14.  The  features  near 
32cm-*  appear  to  be  anti-resonances  associated  with  acoustic  phonons  at  the  SL 
zone  edge,  at  half  the  frequency  of  the  sharp  zone  center  m  *  1  peaks.  The 
bump  at  100  cm-*  is  due  to  second  order  optical  phonon  difference  scattering 
at  wL0(AlAs)-wLA(GaAs). 

Not  all  samples  studied  show  the  sharp  folded  phonons  and  low  background 
that  we  consider  to  be  indicative  of  high  quality  growth.  Figure  15  shows 
data  from  three  samples  whose  Raman  spectra  are  dominated  by  disorder- 
activated  scattering.  The  lower  two  curves  are  from  samples  (14A:0.5,  14A:1) 
and  (38A:0.25,  13A:1),  which  contain  no  pure  GaAs  layers.  DALA  and  DATA 


structure  is  also  seen  in  samples  with  thick  AlAs  layers,  as  in  Fig.  16. 
implies  that  the  disorder  which  relaxes  the  momentum  conservation  does  not 


arise  solely  within  alloy  layers,  but  may  be  more  generally  related  to  the 
difficulty  of  growing  high  quality  layers  with  high  A1  content. 

B.  Optic  Region  (220-420cm~* ) 

Due  to  the  two-mode  behavior  of  AlxGa.xAs  alloys,  the  effect  of  layering 
on  phonons  near  the  GaAs  optical  frequency  in  GaAs-A^Ga^^As  SL's  will  depend 
strongly  upon  the  value  of  x  in  the  barrier  layers.  For  x  <  1  both  layers  may 
support  GaAs-like  optical  vibrations,  resulting  in  a  dispersive  u  vs.  q  with 
small  gaps.  In  this  case,  the  normal  mode  may  have  a  different  number  of 
nodes  in  the  two  layers.  For  AlAs-like  vibrations,  and  for  GaAs-like  vibra¬ 
tions  in  samples  with  x  “  1,  the  phonons  are  highly  damped  in  the  non-propa¬ 
gating  layer.  As  indicated  in  Fig.  2,  these  damped  modes  are  optic-like 
[ Re ( q )  =  x/(a/2)J  in  AlAs  layers  and  acoustic-like  (q  purely  imaginary)  in 
GaAs  layers.  The  localization  in  alternate  layers  produces  flat  optical 
branches.  Whereas  Brillouin  zone  folding  adequately  describes  the  effect  of 
SL  layering  on  the  acoustic  phonons,  the  optic  branches  are  better  regarded  as 
spatially  quantized,  analogous  to  the  quantization  of  electrons  in  potential 
wells. 

Figure  17  shows  data  from  four  GaAs-AlAs  samples  with  GaAs  layer  thickness 
increasing  from  top  to  bottom.  The  dominant  peaks  occur  in  the  (x,y)  geometry, 
indicating  they  have  B2  symmetry.  TO  phonons  occur  in  both  (x,x)  and  (x,y)  at 
~  265cm~l  and  ~  357cm“^.  In  the  upper  two  curves,  weak  Aj  symmetry  peaks  are 
seen  in  (x,x).  These  two  samples  have  been  studied  under  resonant  conditions. 
Figure  18  shows  a  more  detailed  spectrum  of  sample  (20A,  12A:1),  taken  near  the 
band  gap.  Two  peaks  are  seen  in  each  geometry,  in  addition  to  the  TO  phonon. 


ffcsre 


The  resonant  behavior  of  these  phonons  is  given  in  Fig.  19.  The  data  were 
taken  with  =  6765A  by  varying  the  temperature  between  100K  and  300K.  The 
temperature  dependence  of  all  peak  frequencies,  -0.015cm“ */K,  agrees  with  that 
of  GaAs.  The  luminescence  peak  at  an  energy  of  1.805  eV  at  300K  was  used  to 
monitor  the  change  in  band  gap.  Both  the  TO  phonon  and  peak  3  show  little 
resonance  enhancement.  Peaks  2  and  4  resonate  together,  reaching  a  maximum 
20  meV  above  that  of  peak  5. 

Sample  (14A,  12A:1)  shows  a  similar  optical  phonon  structure  and  resonance 
behavior,  except  that  peaks  2  and  4  shift  to  higher  frequencies  with  lower 
exciting  energies  (as  the  laser  is  tuned  through  the  band  gap),  and  their 
resonance  maximum  is  about  50  meV  higher  than  that  of  peak  5. 

A  simple  prediction  of  the  quantized  frequencies  is  difficult  for  lack  of 
a  good  model.  An  upper  limit  for  confinement-induced  shifts  can  be  obtained 
by  regarding  the  bulk  LO  dispersion  as 


2  2  2  2 
oj  =  ui  -  v  q 
L  sn 


where  vg  is  the  longitudinal  sound  velocity.  This  is  the  small  q  limit  of  a 
linear  chain  model,  but  it  overestimates  the  dispersion.  The  folded 
frequencies  are  then  given  by 


222  rAw>2 
U)  *  a)  -  v  (— J 
L  s  vdj' 


where  l  is  the  folding  order  and  dj  is  the  single  layer  thickness.  Using  the 
exact  linear  chain,  we  find  the  result  to  be  somewhat  less  dispersive,  but  it 
is  still  not  accurate  enough,  giving  an  uncertainty  of  several  cm*"*  depending 
on  the  fitting  parameters.  The  existing  neutron  data  on  GaAs  L0  phonon 
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dispersion1* disagree  with  Raman  results  by  6cm-1  at  q  -  0;  they  imply  a 
nearly  flat  LO  branch  out  to  one  third  of  the  bulk  zone  edge. 


I 

Figure  20  is  a  plot  of  observed  optical  phonon  frequencies  in  the  binary 
layer  samples.  The  straight  lines  are  from  Eq.  (52),  using  bulk  values  of  v0 
and  for  several  values  of  £•  Note  that  u>^  is  plotted  vs.  inverse  squared 
layer  thickness.  Peaks  not  otherwise  indicated  appear  in  the  l^-allowed  (x,y) 
geometry.  Peaks  labeled  by  E  have  been  previously  assigned  to  E-symmetry 
phonons,  as  mentioned  below.  Points  along  the  lower  edge  of  both  regions 
corresond  to  TO  phonons. 

The  thickness  dependence  of  the  highest  B£  mode  frequencies  is  rather 
well  described  by  the  i  =  1  line.  These  include  the  peak  labeled  5  In  Fig. 

18.  The  dotted  arrows  suggest  that  the  deviation  shown  by  sample  (llA,  10A:1) 
from  this  line  may  be  due  to  an  incorrect  determination  of  its  relative  layer 
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thicknesses.  Keeping  d  constant  but  adjusting  takes  the  triangles  into 

the  crosses.  X-ray  diffraction  provides  a  precise  determination  of  super- 
lattice  period  from  the  satellite  spacing.  The  measurement  of  individual 
layer  thicknesses  can  be  done  by  Fourier  analysis  of  satellite  intensities  or 
by  estimating  the  average  aluminum  content  from  the  average  lattice  constant. 
The  x-ray  data  are  sometimes  not  good  enough  to  do  this  accurately.  As 
pointed  out  in  connection  with  Fig.  14,  a  good  understanding  of  phonon  folding 
can  allow  Raman  scattering  to  provide  this  structural  Information. 

The  remaining  peaks  between  the  highest  frequency  B2  mode  and  the  TO  mode 
are  difficult  to  assign.  A  rough  correlation  exists  between  the  l  -  2  line 
and  the  second  B2  peak.  Consideration  of  the  mode  symmetries,  and  the  fact 
that  Eq.  (52)  overestimates  the  dispersion,  suggest  that  the  mode  labeled  4 
in  Fig.  18  should  corresponds  to  A  -  2.  The  second  B2  peaks,  labeled  3  in  the 
previous  figure,  then  corresponds  to  £  -  3. 


Previous  reports^  have  suggested  the  E  symmetry  modes  propagating  parallel 
to  the  interfaces  may  appear  in  backscattering  spectra  due  to  disorder.  Such 
modes  have  yet  to  be  systematically  verified,  although  such  scattering  may 
contribute  to  the  structure  typified  by  Fig.  18.  No  additional  optical  phonon 
peaks  were  seen  in  any  of  the  alloy-barrier  samples  included  in  this  study,  as 
were  seen  in  earlier  studies.  This  may  be  due  to  the  high  quality  of  the 
present  interfaces. 

The  interfaces  in  sample  (4lA,  8A:1)  were  intentionally  blurred  by 
annealing,  and  the  optical  frequency  data  are  presented  in  Figs.  21  and  22.  The 
curves  were  taken  in  the  (x,y)  geometry  and  are  normalized  to  50  counts/s  at 
265cm-*  (TO  phonon).  Pieces  of  the  sample  were  annealed  for  0.25,  1,  4,  9  and 
16  hours  at  850°C.  The  L0  phonons  are  seen  to  shift  toward  lower  frequency, 
and  a  low  frequency  tail  appears  on  the  TO  phonon.  The  data  are  approaching 
the  average  alloy  spectrum  with  increasing  annealing.  For  an  average  A1 
concentration  of  x  =  0.16,  confirmed  in  this  sample  by  x-ray  dif fraction, 
the  alloy  LO  peaks  are  expected  at  285  and  370cm-*. *^  The  LO  phonon  in  Fig. 

37  has  moved  below  this  value  to  281cm-*,  perhaps  indicating  a  super lattice 
effect.  Even  after  16  hours  the  composition  profile  should  remain  sine¬ 
like.  22)  No  trace  of  the  behavior  of  an  x  =  .5  alloy  is  seen,  which  might  be 
expected  from  equal  interpentration  at  an  interface.  This  would  give  peaks  at 
271  and  388cm-*. 

In  summary,  we  have  examined  several  models  for  phonon  behavior  in  semi¬ 
conductor  SL's,  and  have  compared  these  with  Raman  scattering  data  from  a 
number  of  aluminum  gallium  arsenide  based  SL  samples.  The  models  describe 
well  the  data  in  the  acoustic  phonon  region.  A  photoelastic  approach  describes 
the  Raman  intensities  from  folded  phonon  doublets.  The  Raman  data  are  shown 
to  provide  structural  information  similar  to  that  gained  from  x-ray 
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diffraction.  In  the  optic  region,  the  correspondence  between  the  models  and 
data  is  not  as  good.  The  behavior  of  the  optic  modes  seems  to  be  best 
described  in  terms  of  quantized  levels,  similar  to  confined  electrons,  due  to 
confinement  within  one  type  of  layer.  This  confinement  depends  of  course  on 
the  alloy  concentration  and  phonon  frequency,  due  to  two-mode  behavior. 

Better  theories,  which  include  the  effects  of  electric  fields,  are  needed  to 
describe  the  optical  behavior. 
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FIGURE  CAPTIONS 

1  Model  used  for  linear  chain  calculations.  AlAs  layer  thickness  is  d2  *  ne . 

2  Real  and  imaginary  part  of  phonon  wavevector  q.  Dashed  lines:  transverse 
modes.  Solid  lines:  longitudinal  modes.  Parameters  of  model  are  fit  to 
circled  points. 

3  Linear  chain  model  calculation  of  phonon  dispersion  curves  for  longitudinal 
and  transverse  modes.  Large  zone:  bulk  GaAs  and  AlAs.  Small  zone:  (5,4) 
superlattice,  q  given  in  A-*. 

4  Frequency  of  interface-like  optic  modes  versus  AlAs  layer  thickness  - 
dj :GaAs  layer,  d2 :AlAs  layer. 

5  Raman  spectrum  of  sample  (42A,  8A:0.3);  A^  “  5145A,  T  *  300K.  Inset  is 
Rytov  model  calculations  showing  q  of  folded  LA  peaks. 

6  Samples  (20A,  1 2A  :  1 )  ,  (4 1 A ,  8A:1),  (59A,  20A:0.3),  and  (4lA,  4lA:0.3); 

*  5145A,  T  »  300K.  Sharp  peaks  are  folded  LA  phonons,  2XA  is  at 
160  cm-*. 

7  Phonon  dispersion  near  first  folded  LA  gap  for  a  (5,4)  superlattice, 
capculated  with  linear  chain  model.  Data  points  are  folded  peaks  in 
sample  (14A,  1 2A : 1 )  at  several  laser  wavelengths.  qz  is  in  units  of  w/d. 

8  Raman  spectrum  of  sample  (14A,  12A:1);  A^  ■  5145A,  T  -  300K,  near  first 
folded  gap.  Strong  peaks  are  LA  phonons,  T  Indicates  transverse  phonons. 


Composite  plot  of  folded  acoustic  phonon  frequencies  vs*  inverse  period. 
Circles  refer  to  samples  with  alloy  barriers,  squares  to  those  with  pure 
AlAs  barriers.  The  two  doublets  at  lower  right  are  transverse  modes. 
Shaded  regions  are  bounded  by  folded  AlAs  and  GaAs  bulk  sound  waves 
(d  -  (2itm/d)vs. 

Raman  spectrum  of  sample  (4lA,  4lA:0.3);  A^  “  5145A,  T  ■  300K.  Inset  is 
Rytov  model  calculation  showing  q  of  folded  LA  phonons  and  missing  second 
order  peaks. 

Raman  spectra  of  sample  (8SA,  88A:0.3)  for  two  values  of  AL;  (x,x) 
geometry,  T  «  300K.  Inset  is  Rytov  model  calculation  showing  q  of  peaks 
in  lower  curve  and  missing  modes.  Upper  resonant  spectrum  shows 
interference-type  behavior. 

X-ray  diffraction  spectrum  of  sample  (85A,  88A:0.3)  showing  negligible 
second  order  satellite. 

Raman  spectra  of  folded  LA  phonons  from  sample  (4lA,  8A:1)  annealed  at 
850K  for  0,  1,4,  and  9  hrs;  AL  ■  5145A,  T  -  300K,  (x,x)  geometry. 

Raman  spectra  near  the  band  gap  of  sample  (14A,  12A:);  AL  -  6328A ,  T  ■ 
300K.  Spikes  above  100  cm“*  are  Ne  calibration  lines. 

Raman  spectra  of  3  samples  showing  disorder  activated  structure.  From  top 
to  bottom  the  samples  are  (llA,  10A:1),  (38A:0.25,  13A:1),  and  (14A:0.5, 

1 4 A  :1)  A.  -  5145A,  T  -  300K. 


16  Raman  spectra  of  sample  (12. 5A,  37.5A:1)  showing  disorder  activated 
structure  at  two  temperatures;  X^  =  5145A. 

17  Optical  phonon  Raman  spectra  of  4  samples  with  AlAs  barriers.  Strong  peaks 
are  in  (x,y)  geometry,  weak  curve  is  (x,x).  From  top  to  bottom  the  samples 
are:  (14A,  12A:1);  (20A,  12A:1);  (27A,  llA:l);  (4lA,  8A:1). 

18  Raman  spectrum  near  GaAs-like  optic  branch  of  sample  (20A,  1 2A  :  1 ) .  X^  * 
6764A,  T  =  180K. 

19  Resonance  behavior  of  peaks  labeled  in  Fig.  18.  A  refers  to  lowest  folded 
acoustic  phonon.  Band  gap  was  tuned  through  X^  by  varying  temperature. 
Intensity  not  corrected  for  absorption.  uQ  is  300K  luminescence  peak. 

20  Composite  plot  of  optical  phonon  frequencies  squared  vs.  inverse  square 

layer  period,  d^  =  GaAs  layer  thickness,  ^  “  AlAs.  Straight  lines  are  u>^ 

7  7  7 

=  -Vg(Ax/d^  2)  for  A  *  1,2,  or  3.  Unless  specified,  peaks  appear  in 
(x,y)  geometry. 

21  Raman  spectra  of  GaAs-like  optic  phonons  from  sample  (4lA,  8A:1)  annealed 
at  850K  for  0,  1,4,  16  hours;  XL  -  5145A,  T  «  300K,  (x,y)  geometry. 

22  Same  as  Fig.  21  for  AlAs-like  optic  region. 
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The  first  reported  Raman  scattering  experiments  have  been  performed  on 
single  crystal  and  amorphous  films  of  the  metastable  alloy  (GaSb)^  x^ex 

with  compositions  across  the  pseudobinary  phase  diagram.  In  crystalline 
films,  the  optical  phonons  exhlhlt  a  "one-two"  type  mode  behavior. 

Broadening  and  softening  of  the  Raman  peaks  with  Increasing  alloy  con¬ 
centrations  are  attributed  to  a  relaxation  of  q-vector  selection  rules 
due  to  substitutional  disorder  on  both  the  cation  and  anion  sublattices. 

Additionally,  disorder  Induced  scattering  from  the  zone-edge  acoustic 
phonons  was  observed.  In  amorphous  alloy  films,  the  reduced  Raman 
spectra  were  compared  to  the  one-phonon  densities  of  states  of  the  end- 
member  crystals.  The  .result Ing  apparent  lack  of  polarization  dependence 
Indicated  that  the  amorphous  films  were  of  the  random  network  type. 


Semiconducting  alloys  have  been  extensively 
studied,  not  only  because  of  their  practical  de¬ 
vice  applications, ‘but  because  of  fundamental 
Interest  in  disordered  systems.  Tn  particular, 
the  lattice  dynamics  of  mixed  crystal  and  amor¬ 
phous  systems  have  been  of  great  interest.  Raman 
and  Infrared  spectroscopies  have  been,  employed 
to  Investigate  the  behavior  of  long-v.-avelength 
optical  phonons  in  nixed  crystals  as  a  function 
of  composition.  The  observed  behavior  generally 
falls  Into  one  of  two  classes. The  spectra  of 
mixed  crystals  with  one-mode  behavior  display  a 
single  set  of  optical  modes;  as  the  composition 
Is  varied  the  strength  of  these  modes  remains 
approximately  constant,  while  their  energy  shifts 
continuously  from  the  optical  mode  energy  of  one 
end-member  component  to  that  of  the  other  com¬ 
ponent.  This  pattern  is  typically  found  In  I -VI I 
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materials.  Mixed  crystals  displaying  a  two-mode 
behalvor  possess  two  sets  of  optical  modes  oc-' 
currlng  close  to  those  of  each  component  with  the 
strength  of  each  mode  dependent  on  the  concentra¬ 
tion  of  the  corresponding  component.  Many  of  the 
III-V  ternary  alloys  exhibit  this  two-mode  £ehav- 
lor.  However,  several  semiconductor  alloys  pres¬ 
ent  more  complex  patterns.  Ternary  alloys  such 
as  InA8i_xSbjj  and  Ga^_xInxSb-^  are  found  to  dis¬ 
play  an  Intermediate  behavior  with  two  sets  of 
optic  modes  observed  over  pert  of  the  composition 
range,  with  only  one  mode  seen  over  the  remaining 
range  of  compositions  ("one-two"  mode  behavior). 
The  alloy  system  Slj  Ge  exhibits  a  three  mode 
behavior;4'3  because  Botn  constituents  can  sub¬ 
stitute  freely  on  the  diamond  lattice,  Sl-Sl,  Ge- 
Ge,  and  Sl-Ge  bonds  are  present,  and  the  modes 
observed  are  attributed  to  this.  Several  theo¬ 
ries  have  been  advanced  to  predict  and  explain 
the  long-wavelength  mode  behavior  In  mixed  crys-* 
tale.6-" 

In  this  paper  we  report  the  first  Raman  re¬ 
sults  on  a  novel  semiconducting  alloy  (GaSb)x-x 


896 


Vol.  45,  No.  10 


CRYSTALLINE  AND  AMORPHOUS  (GaSb) ,_xGex  SEMICONDUCTING  FILMS 


Ge*  which  has  recently  been  grown  in  both  cry¬ 
stalline10*11  and  amorphous1'  forms  with  concen¬ 
trations  ranging  across  the  pseudobinary  phase 
diagram.  It  should  be  noted  that  this  mixed 
crystal  differs  in  several  ways  from  previously 
studied  alloys.  The  single  crystal  alloys  are 
metastable  over  most  of  the  Investigated  concen¬ 
tration  range.  The  pseudobinary  phase  diagram 
of  this  system13  shows  that  the  equilibrium  so¬ 
lid  solubilities  of  GaSb  In  Ge  and  Ge  in  GaSb 
are  only  6.5  ±  0.5  and  2.0  ±  0.5  mole  X,  re¬ 
spectively.  Another  difference  from  ternary 
alloys  Is  that  substltutuional  disorder  exists 
on  both  sublattices  of  th'e  zincblende  or  diamond 
structure.  It  Is  unique  that  the  substituting 
species  are  not  isovalent  substitutions;  Ge  pro¬ 
vides  both  acceptors  and  donors  when  it  is  an 
impurity  in  GaSb,  since  it  occupies  both  cation 
and  anion  sites  in  the  lattice,  while  GaSb  im¬ 
purities  in  Ge  provide  Ga  acceptors  and  Sb 
donors. 

All  (GaSb^.xGe^  alloys  were  deposited 
using  a  multitarget  radio-frequency  sputtering 
system  that  has  been  described  elsewhere. 14’ 10 
The  single  crystal  films  were  epitaxially 
grown  on  semi-insulating  (100)  GaAs  substrates 
under  excess  antimony  pressures  at  growth  tem¬ 
peratures  ranging  from  470  to  485°C.  The 
amorphous  (GaSb)i_xGex  mixtures  were  grown  on 
Corning  7059  glass  substrates  at  60°C  with  no 
antimony  overpressures.  In  both  cases,  the 
average  film  thickness  was  1.5  pm. 

Film  compositions  were  determined  from  the 
wavelength  dispersive  analyses  in  a  JEOL  elec¬ 
tron  microprobe  with  an  accuracy  of  0.5  at.Z. 

The  structural  perfection  of  the  crystalline 
films  was  confirmed  with  the  use  of  x-ray  dif- 
fractometry,  Debye-Scherrer  patterns,  and  elec¬ 
tron  channelling.  X-ray  diffraction  spectra 
exhibited  very  sharp  (400)  peaks  whose  widths 
were  limited  by  instrumental  resolution. 

Electron  channelling  patterns,  taken  with  a  10 
um  aperture,  were  sharp  and  found  to  be  invari¬ 
ant  while  scanning  over  the  sample,  demonstrat¬ 
ing  the  homogeneity  of  the  film.  The  film 
lattice  constants  obeyed  Vegard's  law,  varying 
linearly  with  composition. 10  Post -annealing 
experiments  showed  that  the  single  crystal  meta¬ 
stable  films  exhibited  good  temporal  and  thermal 
stability.  The  transformation  to  the  equilibri¬ 
um  two  phase  state  was  found  to  occur  at  temper¬ 
atures  ranging  from  475  to  525°C,  depending  upon 
the  alloy  composition.  From  measurements  of  the 
exothermic  hept  of  transformation  as  well  as  cal¬ 
culations  of  the  total  free  energy  change,  the 
lifetime  of  the  single  phase  metastable  state  at 
room  temperature  was  estimated  to  be  of  the  order 
of  1029  years.  In  the  present  Raman  experiments, 
irradiation  with  a  well  focused  one  watt  Ar+ 
laser  beam  had  no  effect  on  sample  integrity. 

The  Raman  measurements  were  performed  with 
the  sample  at  room  temperature  and  under  vacuum 
to  avoid  spurious  scattering  from  the  air  mole¬ 
cules.  Light  at  a  wavelength  of  514.5  nm  from  an 


Ar  laser  was  incident  perpendicular  to  the 
(100)  surface  of  the  sample  with  the  scattered 
light  collected  in  the  true  backscatterlng  geome¬ 
try.  A  double  monochromator  and  conventional 
photon  counting  electronics  were  used.  The  laser 
power  was  one  wstt  for  the  single  crystal  alloy 
films.  However,  in  the  amorphous  films  the  power 
was  limited  to  20  mW  in  order  to  ovoid  spontane¬ 
ous  crystallisation. 17 


Figure  1  shows  the  Raman  spectra  from  six 
single  alloy  samples,  pure  Ge,  and  pure  GaSb. 
These  spectra  were  taken  with  the  Incident  light 
polarized  along  the  [010]  crystal  axis  and  the 
scattered  light  analyzed  along  the  [001]  axis. 


Fig.  1  Raman  spectra  from  single  crystals: 
a)  Ge  b)  (GaSb)0  27Ce0  „ 

c)  (GaSb)0  A4Ge0  56  d)  (GaSb>0  ^Ge^ 
e)  (GaSb)0  76Ce0  24  f)  (Ca3b)Q >86Ce0 
g)  (GaSb)0  94Ge0  Q6  h)  GaSb.  All 

spectra  were  taken  in  the  z(xv)z  con¬ 
figuration  with  a  514.3  nm  Ar*  laser 
beam  at  room  temperature. 

According  to  the  usual  Raman  selection  rules,18 
in  this  configuration  the  GaSb  L0  mode  and  the  Ge 
optic  mode  are  allowed,  but  the  GaSb  TO  mode  is 
forbidden.  Figure  1  shows  that  two  optic  modes 
were  observed  in  the  alloy  spectra,  one  near  the 
GaSb  L0  mode  at  234  cm-1  and  one  near  the  Ge  op¬ 
tic  mode  at  300  cm-1.  The  GaSb-llke  mode  did 
not  shift  appreciably  in  frequency  as  the  GaSb 
content  was  decreased;  hovever,  its  intensity 
fell  rapidly  and  the  peak  broadened  considerably 
so  that  it  was  just  discernible  la  the  (GeSb>o.27 
Geo. 73  spectra,  (hi  the  other  hand,  the  frequency 
of  the  Ge-llke  node  decreased  from  300  to  257 
cm-1  and  the  Intensity  decreased  in  an  approxi¬ 
mately  linear  fashion  as  the  Ge  concentration 
was  dacreased  from  100  to  6  mole  X.  The  line- 
width  of  this  node  remained  almost  constant  in 
the  alloys. 
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Figure  2  shows  the  frequency  of  both  the 
GaSb-like  and  Ge-llke  modes  versus  alloy  concen¬ 
tration.  The  strong  damping  of  the  GaSb-llke 
mode  with  decreasing  GaSb  concentration  suggests 
that  (GaSb)j_xGex  follows  an  Intermediate  one- 


Fig.  2  Long-wavelength  optical  phonons  in 

(CaSb)^_xGgx  versus  alloy  concentration. 

two  mode  pattern.  In  two-mode  behavior,  the  op¬ 
tical  modes  of  the  pure  crystal  jend-members 
evolve  Into  a  localized  impurity  .mode  of  that 
constituent  in  the  other  end-member.  Intermedi¬ 
ate  behavior  occurs  when  creation  of  a  distinct 
localized  mode  does  not  take  place.  In  this  al¬ 
loy,  we  suggest  that  small  concentrations  of  Ga 
and  Sb  in  Ge  do  not  induce  localized  modes  In 
the  lattice  for  the  following  reasons:  (1)  Ga 
has  too  small  a  mass  difference  with  Ge  to  sig¬ 
nificantly  perturb  the  lattice,  and  (2)  the  Ge 
crystal  does  not  have  a  gap  in  Its  phonon  densi¬ 
ty  of  states  for  creation  of  the  gap  mode  associ¬ 
ated  with  heavy  substitutional  Sb  Impurities 
(Msb^fce  “  1.69)'.  Spectra  taken  from  a  99.8 
mole  Z  Ge  sample  showed  no  evidence  for  Ga  or  Sb 
related  local  modes.  In  contrast,  at  the  other 
end  of  the  alloy  concentration  range,  the  Ge-. 
like  mode  merged  into  a  local  modeat  257  cm-1 
as  shown  in  Figure  1(g).  We  attribute  this  to 
Ge  atoms  residing  on  Sb  sites.  Ge  atoms  on  Ga 
sites  have  too  small  a  mass  difference  for  cre¬ 
ation  of  a  gap  mode. 

The  strong  broadening  of  the  optical  modes 
and  the  appearance  of  low-frequency  modes  in  Fig. 
1  are  attributed  to  first-order  scattering  in¬ 
duced  by  the  substitutional  disorder  in  the  mixed 
crystals.  This  disorder  relaxes  the  qM)  Raman 
selection  rule  and  allows  phonons  away  from  the 
zone  center  to  become  Raman  active.  The  broad 
feature  found  near  150  cm-1  corresponds  to  the 
zone-edge  longitudinal  acoustic  phonon  frequen¬ 
cies  in  GaSb  and  Ge;  structure  seen  at  about  50 
cm*1  matches  the  zone-edge  transverse  acoustic 
phonon  frequencies.  These  observed  peaks  are  in¬ 
terpreted  as  disorder  activated  longitudinal 
acoustic  and  transverse  acoustic  modes.  The  re¬ 


laxation  of  the  wave  vector  selection  rule  also 
affects  the  width  and  frequency  of  the  optical 
phonon  peaks.  Off-zone  center  optical  phonons 
also  become  Raman  active  due  to  the  disorder. 
This  broadens  the  peak,  down-shifts  the  peak  po¬ 
sition  in  energy,  and  gives  the  peak  an  assyme- 
trlc  appearance,  because  of  the  negative  slopes 
in  the  optical  phonon  dispersion  curves  of  both 
Ge  and  GaSb.  In  spite  of  the  substitutional 
disorder  In  these  crystalline  samples,  the  po¬ 
larization  dependence  of  the  optical  modes  was 
the  same  as  that  of  the  zone-center  modes  ob¬ 
served  in  pure  GaSb  and  Ge.  The  disorder-al¬ 
lowed  acoustic  phonons  showed  little  dependence 
on  the  polarization.  This  is  expected,  since 
the  frequencies  of  phonons  from  points  such  as 
X,  L,  and  W  are  nearly  the  same,  and  a  superpo¬ 
sition  of  spectra  would  contain  all  Raman-active 
symmetries. 

Spectra  obtained  from  amorphous  (GaSb) j.x 
Gey  films  are  shown  in  Fig.  3.  In  amorphous  ma¬ 
terials  the  q-vector  selection  rule  breaks  down 


Fig.  3  Raman  spectra  of  (GaSb)1_xGex  amorphous 

alloys  and  both  end  members:  a)  GaSb 
DOS,  from  Ref.  A,  convoluted  with  a 
Gaussian  FWHM  ”  20  cm“l 

b)  a-OaSb 

c)  (CaSb)0  90C.e0  10 

d)  (CaSb). >56Gc0 

e>  (CaSb)0<13Ce0  8? 

f)  a-Ce 

g)  Ge  DOS,  from  Ref.  B,  convoluted awith 
a  Gaussian  FWH.M  *  50  cm-* 

A)  Marvin  K.  Farr,  Joseph  G.  Taylor  and  S.  K. 
Sinha,  Physical  Review  BU,  1587  (1973). 

B)  G.  Nelin  and  G.  Klllson,  Phvsical  Review 
B5  3151  (1972). 
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completely,  and  the  Raaan  spectra  ara  a  measure 
of  the  phonon  density  of  aCatea  (DOS).  An  appro¬ 
ximate  expresalon  for  the  Raaan  intensity  I(u)  In 
thla  case  la:19 


E  Cb(»)pb(a) 


n(m.T)  +  1 


where  n(u,T)  la  the  Boae-Elnatein  distribution 
function,  p.  (ai)  Is  the  vibrational  density  of 
states  In  band  b,  and  C.  (u)  la  the  coupling 
strength  between  the  exciting  light  and  the  phonon 
band.  The  data  in  Fig.  3  are  plotted  as  redu¬ 
cing  Raaan  spectra: 

Vw)  *  n(a,T>  +  lI(w)  “p(“> 

where  the  proportionality  to  the  total  density 
of  states  holds  If  we  assuae  that  Cj,(m)  la  con¬ 
stant  with  u>  and  b.  For  comparison  with  Ij^w), 
broadened  phonon  densities  of  states  for 
GaSb  and  Ge  are  also  plotted  in  Fig.  3.  These 
plots  were  generated  by  convoluting  a  Gaussian 
with  the  one-phonon  density  of  states,  choosing 
the  half  width  of  the  Gaussian  equal  to  the  en¬ 
ergy  difference  between  the  LO  phonon  and  the 
cutoff  of  first  order  scattering  In  the  amor¬ 
phous  spectra.18  For  alloy  samples  with 
x  <  .87,  the  spectra  contain  contributions  from 
the  DOS  of  both  end  members.  The  sample 
(GaSb) _ 1 3Ge<  gy  (Fig.  3e)  does  not  display  any 
contribution  from  GaSb,  consistent  with  the 


one-two  mode  pattern  postulated.  The  difference 
between  the  spectra  and  the  DOS  of  the  pure  com¬ 
pound  a,  particularly  In  the  Intensities  of  fea¬ 
tures  at  low  u,  la  moat  probably  due  to  vari¬ 
ations  in  (^(w) .  There  was  no  apparent  polar¬ 
ization  dependence  of  the  spectra,  leading  to 
the  conclusion  that  the  amorphous  f 11ms  have  a 
continuous  random  network-type  structure.21 

In  summary,  Raaan  scattering  results  on 
(GaSb)  j-xGCx  single  crystals  show  an  Intermedi¬ 
ate  one-two  optic  phonon  mode  behavior.  This  is 
attributed  to  the  lack  of  a  localized  mode  by 
GaSb  Impurities  In  Ge.  A  local  mode  is  observed 
In  GaSb  due  to  Ge  Impurities  on  Sb  sites  in  the 
GaSb  lattice.  Substitutional  disorder  In  the 
crystalline  alloy  Induced  substantial  broadening 
and  softening  of  optical  modes  In  addition  to 
disorder -Induced  scattering  from  the  zone-edge 
acoustic  phonons.  It  is,  however.  Interesting 
to  note  that  the  relaxed  Raman  selection  rules 
persist  In  the  crystalline  alloy  films,  whereas 
no  polarization  dependence  was  found  in  the  amor¬ 
phous  forms. 
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A  great  deal  of  effort  has  gone  into  the  design  and  construction  of  a 
tandem,  multipass  Fabry- Perot  interferometer  system  for  use  in  studying 
bulk  and  surface  acoustical  phonons  in  semiconductor  alloys  and  super¬ 
lattices.  The  system  is  still  not  operational;  work  on  it  will  continue 
using  other  research  support.  Improvements  to  an  old  far-infrared  Fourier 
interferometer  have  been  made,  and  the  instrument  will  be  used  to  study 
optical  phonons  in  semiconductor  alloys  and  superlattices  in  future  work. 

Raman  studies  have  been  made  on  semiconducor  alloys  and  superlattices. 
One  reprint  and  one  preprint  of  papers  on  (GaSb).  Ge  metastable  alloys 

are  included.  A  great  deal  of  data  has  been  taken  on  the  GaSb.  As 

l-x  x 

metastable  alloy  system,  but  analysis  is  not  yet  complete.  An  extensive 
study  has  been  completed  on  folded  acoustic  and  quantized  optic  phonons 
in  GaAs/Ga1_xAlxAs  superlattices,  and  a  paper  has  recently  been  submitted 
to  Physical  Review  B.  A  preprint  is  enclosed. 
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Preprint  of  paper  to  be  presented  at  the 
17th  International  Conference  on  the  Physics  of  Semiconductors 
San  Francisco,  California,  August  6-10,  1984 
Proceedings  to  be  published  by  Springer-Verlag 
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In  an  extension  of  earlier  work  (1)  we  have  used  Raman  Scattering 
to  reveal  structural  information  in  mixed,  metastable,  crystalline 
[GaSb] i-xGe2x”lGai-xGex] [Sbi_xGex] .  These  alloys  were  deposited  on 
GaAs  substrates  using  a  multitarget  r.f.  sputtering  system  (2). 

Room  temperature  measurements  were  performed  under  vacuum,  with 
the  incoming  laser  beam  perpendicular  to  a  <100>  plane,  and  photon 
polarizations  along  (100]  directions.  Figure  1  shows  the  Raman  spectra 


Fig.  1:  Room  temperature,  Raman 
scattering  of  mixed 
(GaSb) i-xGe2x*  Scattering 
geometry:  z(x,y)z. 


of  some  of  the  mixed  crystals  and  of 
the  pure  GaSb  and  Ge  components  (all 
normalized  to  same  maximum  peak 
height) .  When  small  concentrations 
of  GaSb  are  introduced  into  Ge  or 
small  conentrations  of  Ge  introduced 
into  GaSb,  the  zone-center  phonon 
frequency  shifts,  broadens,  and 
becomes  asymmetric.  Figure  2  plots 
the  frequencies  versus  "percent 
Germanium”,  which  we  define  as  lOOx. 
For  x<0.1  the  GaSb-like  phcnon  shows 
a  TO-LO  splitting,  which  cannot  be 
detected  for  x>0.1.  The  frequency 
of  the  Ge-like  optical  phonon 
increases  almost  linearly  with  x, 
whereas  the  frequency  of  the 
GaSb-like  phonon  remains  almost 
constant  for  0.3  <  x  <  0.75.  For 
x>0.75  this  mode  can  no-longer  be 
detected,  consistent  with 
expectations  for  behavior  of 
isolated  substitutional  Ga  and  Sb 
atoms  in  Ge:  In  low  concentrations 
neither  atom  will  give  a  local  mode 
or  a  gap  mode.  Ge,  however,  does 
produce  a  local  mode  as  a  dilute 
impurity  in  GaSb. 
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Pig.  2:  Frequency  dependence  of 
peaks  in  (GaSb)i_xGe2X  as  a 
function  of  Ge  concentration. 
Error  oars  are  qualitative 
estimates. 
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Fig.  3:  Width  and  asymmetry  at 
half  maximum  of  the  optical 
phonons . 


The  half  widths  of  the  optical 
phonons  (at  half-maximum)  will  be 
denoted  by  ra  and  below  and  above 
the  peak,  respectively.  The  total 
width  r  =  ra+rt),  and  the  asymmetry 
fa/rb  are  plotted  versus  x  in  Fig. 

3.  For  the  GaSb-like  mode  r  is  an 
increasing  function  of  x  until 
x=0.75,  where  the  mode  can  no  longer 
be  detected,  r  for  the  Ge-like  mode 
rises  first  linearly  with  (1-x) , 
reaches  a  maximum  for  x*0.35  and 
then  falls.  The  maximum  values  of  r 
(GaSb)  and  r(Ge),  47  and  35  cm~l 
respectively,  are  close  to  those  of 
the  pure  amorphous  components 
(S2cm_l  for  a-GaSb  and  47cm- 1  for 
a-Ge).  The  asymmetry  ra/rb  of  the 
GaSb-like  mode  increases  with  x 
until  x*0.25.  For  x>0.25  ra/rb 
decreases,  whereas  r  continues  to 
increase.  Similarly,  the  asymmetry 
of  the  Ge-like  mode  increases  with 
(1-x)  until  (1-x)*  0.20.  For 
(1-x) >0.25  the  asymmetry  decreases, 
whereas  r  continues  to  increase, 
until  (1-x) *0.65.  This  behavior  is 
different  from  that  found  in  the 
quasi-binary  III-V  alloys  Ga}-xAlxAs 
and  Gai_xInxAs,  where  ra/Tb  is  an 
increasing  function  of  r  (3)  and  in 
GaAs  implanted  with  As+  where  ra/rb 
and  r  in  samples  that  remain 
crystalline,  are  increasing 
functions  of  fluence  (4).  In  these 
three  systems  r  was  never  greater 
than  about  14  cm"*,  and  the  maximum 
value  of  ra/rb  was  about  2.2. 

By  its  very  nature  [Gai_xGexl 
(Sb2-xGexl  has  more  disorder  than, 
say,  GaSbi-xAsx,  since  in  the 


present  case  both  sublattices  of  the  underlying  lattice  become 
disordered.  The  greater  values  of  r  observed  here  are  at  least  partly 
due  to  this  effect.  An  important  issue  is  whether  the  non-monotonic 
relation  of  ra/rfc  to  r  is  simply  due  to  the  larger  amount  of  disorder, 
rather  than  to  qualitative  changes  in  the  nature  of  the  disorder.  The 
abrupt  changes  in  ra/Tb  (GaSb)  near  x=0.20  and  ra/rfc(Ge)  near  x=0.80 
seem  to  occur  over  too  small  a  change  in  x  to  be  simply  a  quantitative 
effect  of  disorder. 

In  references  (3)  and  (4)  the  results  were  analyzed  using  the 
"Spatial  correlation  model*  which  intoduces  a  correlation  length  L  to 
account  for  the  wave-vector  relaxation  (5).  The  contribution  of 
phonons  of  finite  q  with  frequency  m(q)  and  width  rQ  is  accounted  for 
by  use  of  the  following  expression  for  the  Raman  line-shape: 


!(<■>) 


3 

_ d_q _ 

[w-iD(q)  ]  2  +  [  rQ/2]  2 
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u(q)  was  taken  to  be  that  given  by  a  linear  chain  model  (3),  (4).  The 
resulting  values  of  ra,  r^,  and  ra/rb  are  increasing  functions  of 
L"1.  A  rough  extrapolation  of  these  results  to  large  r,  suggests  that 
the  T-values  in  Fig.  3  imply  minimum  values  of  L  of  order  10 A.  When  a 
more  realistic  model  is  used  for  u(q),  and  when  both  TO  and  LO  modes 
are  included,  it  is  possible  that  for  very  large  L_1  -  107  cm-1  Eq.  (1) 
will  give  a  line-shape  with  reduced  asymmetry,  but  the  change  will  not 
be  an  abrupt  function  of  L”*.  The  ra/Tb  data  (Fig.  3)  suggest 
qualitative  changes  in  the  disorder  at  x«0.20-0.30  and  x-0.70-0.80. 

In  (GaAs)i-xGe2X,  there  is  experimental  and  theoretical  evidence  for  a 
transition  from  an  average  zincblende  structure  to  an  average  diamond 
structure  as  x  increases  past  xc=>  0.35  (6,7).  X-ray  diffraction  data 
on  (GaSb) 2.xGe2x  show  a  similar  transition  at  xc«0.30  (8).  For  x>xc 
there  is  no  long-range  zincblende  order.  The  Raman  results  do  not 
support  the  picture  that  site  occupancy  is  completely  random  when  x>xc, 
because  if  that  were  the  case  we  should  see  Raman  scattering  from  Sb-Sb 
bonds.  If  these  bonds  were  metallic,  peaks  should  be  seen  at  115  cm“l 
and  150  cm-l  (9).  If  the  bonds  were  covalent,  a  peak  would  be  seen  at 
about  193  cm“l  (a  frequency  slightly  scaled  down  from  the  Raman 
frequency  in  gray  tin  (10)].  There  are  no  sharp  features  near  115,  150 
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or  190  cm“l  that  could  be  identified  with  Sb-Sb  bonds.  This  fact 
together  with  the  observation  that  the  GaSb-like  phonon  continues  to 
exist  up  to  x=0.75  implies  that  Sb-Ga  (and/or  Sb-Ge)  bonds  exist,  but 
not  Sb-Sb  bonds.  If  Sb-Ga  bonds  tend  to  be  favored  over  Sb-Ge  bonds, 
due  to  the  coulomb  interaction,  there  must  be  small  regions  of 
zincblende  structure  alternating  with  regions  of  anti-zincblende 
structure.  The  large  values  of  L-*  suggested  by  the  large  values  of  r 
imply  that  these  regions  would  be  as  small  as  10  A  in  extent.  This 
figure  is  consistent  with  our  failure  to  see  any  microstructure  using  a 
scanning  transmission  electron  microscope  with  a  resolution  of  50  A. 

Whatever  the  nature  of  the  short-range  order  or  microstucture,  the 
ra/rb  data  suggest  that  the  structure  changes  at  x-0.20-0.30  and  again 
at  x=0 .70-0 .80 . 
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